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ABSTRACT 

The statistical mechanical investigation of Violent Relaxation of phase space ele- 
ments of different densities first derived by Lynden-Bell (1967) is re-examined. It is 
found that the mass independence of the equations of motion of Violent Relaxation 
calls for a constraint on the volume of the phase space elements used to formulate the 
statistical mechanical description of Violent Relaxation. In agreement with observa- 
tions of astrophysical objects believed to have been subject to Violent Relaxation (e.g. 
clusters of galaxies), the coarse grained phase space distribution / of the final state 
in the non-degenerate limit turns into a superposition of Maxwellians of a common 
velocity dispersion. Thus, the velocity dispersion problem present in the investigation 
of Lynden-Bell (1967) is removed. 

Subject headings: galaxies: formation — galaxies: interactions — gravitation — meth- 
ods: statistical 
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1. INTRODUCTION 

Since the classical paper by Lynden-Bell (1967), 
the statistical investigation of Violent Relaxation 
has been a topic of crucial interest in the dis- 
cussion of the formation of relaxed astrophysi- 
cal objects like star clusters, galaxies and clus- 
ters of galaxies. Its reformulation by Shu (1978) 
in terms of particles instead of phase space el- 
ements represents a change in terminology and, 
to some extent, a re-interpretation of the Liou- 
ville theorem. However, the basic result of both 
attempts is essentially the same and has been a 
subject of continuous discussion. The debate of 
Madsen (1987) and Shu (1987) concerning the 
foundations of a statistical description of Violent 
Relaxation provides a summary of resolved and 
unresolved issues for both the particle and the 
phase-space representation. 

Starting with Tremaine, Henon, & Lynden- 
Bell (1986), the form and evolution of the H- 
function during Violent Relaxation became a 
subject of growing interest. Ever since one of the 
main lines along which the statistical mechanical 
investigation of violent relaxation proceeds, deals 
with this problem (e.g., Stiavelli & Bertin 1987; 
Kandrup 1987). Only recently (Soker 1996) it 
became clear, that the H- function of collision- 
less self-gravitating systems does not necessarily 
increase monotonically during the time evolution 
of the system. 

In the present paper we deal with another is- 
sue central for both the particle and the phase- 
space description of Violent Relaxation. Since 
the very first investigation by Lynden-Bell (1967), 
it is considered as one of the most important 
flaws of the theory (Madsen 1987; Shu 1987) 
that despite the collisionless nature of Violent 
Relaxation the statistical mechanical approach 
predicts a thermalized final state, i.e. a state 
in which the square of the velocity dispersion of 
phase space elements of different densities is in- 
versely proportional to the density. In order to 
cure this problem, Shu (1987) applied very strin- 



gent assumptions on the initial phase space den- 
sity to his particle approach. 

Here we re-examine the statistical mechanics 
of Violent Relaxation in terms of phase space el- 
ements of different densities. Since the basic ob- 
ject on which the statistical mechanical approach 
of Lynden-Bell (1967) is based are phase space el- 
ements of constant volume but of different mass, 
it is found that the mass independence of the 
equations of motion of Violent Relaxation is not 
fully represented. Considering phase space ele- 
ments of different volume but of constant mass 
removes the velocity dispersion problem tenden- 
tially present in the investigation of Lynden-Bell 
(1967). 

2. VIOLENT RELAXATION OF PHASE 
SPACE ELEMENTS OF DIFFERENT 
DENSITIES 

The collisionless nature of the process of Vi- 
olent Relaxation allows to consider its dynamics 
in terms of the one-particle phase space [i. The 
Vlasov equation 



D t f =<9 x /-v-V x $-cV = 



(1) 



governs the collisionless dynamics of the fine- 
grained phase space density / during Violent 
Relaxation. Violent Relaxation is essentially a 
self-gravitating process. Accordingly, the gravi- 
tational potential $ in (|l|) is given by 



x 



x' 



(2) 



Equations (|T|) and (|2|) describe the dynamics of 
Violent Relaxation in terms of the continuum 
limit without reference to the individual mass of 
the objects involved. We will restrict ourself to 
the continuum limit description and assume its 
validity hereafter. As a consequence, the discus- 
sion is limited to the phase space element ap- 
proach (Lynden-Bell 1967) and does not refer 
to its counterpart, the particle approach of Shu 
(1978) (see also Madsen (1987) and Shu (1987)). 
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By considering Violent Relaxation in the ter- 
minology of phase space elements, the attempt of 
Lynden-Bell (1967) incorporates the constraint 
on the time evolution of the /i-space provided by 
the Vlasov equation (jl]) stating the constancy of 
the phase space density / along trajectories in 
/i-space. As a consequence of the Vlasov equa- 
tion, initially not overlapping phase elements do 
never overlap. Therefore, a microscopic exclu- 
sion principle for phase space elements in //-space 
is established. The second characteristic feature 
of Violent Relaxation reflected by equation (|l|) 
is the independence of the equations of motion 
on the individual phase space element or parti- 
cle masses. It is incorporated into the statistical 
mechanical picture by the dimensionality of the 
/i-space which is chosen to yield the mass density 
in units mass per spatial and velocity volumes, 
A 3 x and A s v, respectively. 

In order to apply statistical mechanical meth- 
ods, the /i-space is divided into a large number 
of microcells. With r\ the phase space density of 
a phase space element (i.e. an occupied micro- 
cell) and co its volume, the mass associated with 
the phase element is rjco. A microstate may then 
be described by the set of occupied microcells Uj 
and their density rjj . While the volume occupied 
by the microstate in the //-space is 

&H = T, N M> ( 3 ) 
j 

the corresponding volume Ar in T-space is 

Ar = n-f J (w) 3Ar ' (4) 

j 

where Nj is the total number of phase space ele- 
ments of density rjj occupying microcells of vol- 
ume ujj. The first crucial point is that the phase 
space volume in T-space (|4|) not only depend on 
the volumes Uj but does also depends on the 
phase space density rjj. The second point con- 
cerns the total energy of the microstate 

^Ew^ + ^l (5) 

j 



where Vj is the velocity associated with the ith 
microceil, and <&j is the gravitational potential 
which the microceil is exposed to. What is of 
importance here is that the total energy E is 
related to the T-space, i.e. it is not possible 
to derive E independently of the absolute value 
of the phase space density r/j, or strictly speak- 
ing, of the mass rjjiVj of the phase space element. 
This dependence is recovered in the macroscopic 
constraints the system is subject to in the sta- 
tistical mechanical approach. In this sense, the 
basic objects on which the Lynden-Bell (1967) 
statistic is based are not the (mass independent) 
phase space element volumes Uj in /i-space but 
their mass dependent counterparts ujj(rjjU!j) 3 in 
T-space. Therefore the discussion does, at least 
with respect to the energy constraint, not reflect 
the universal mass independence of the equations 
of motion of Violent Relaxation. It is indeed easy 
to see that it is the energy constraint which for- 
mally introduces the dependence of the coarse 
grained phase space distribution on the phase 
space density rjj in the approach of Lynden-Bell 
(1967). This means that in order to incorpo- 
rate or represent the mass independence it is not 
enough to just chose the dimensionality of the 
//-space as discussed above. 

We seek for a generalized constraint related 
to the division of the /i-space into phase space 
elements and microcells, respectively. This di- 
vision of the /i-space refers to the intermediate 
step of the determination of the most probable 
state in terms of occupation numbers of macro- 
cells. While this step relies on a specific division 
of the /i-space, the final result is, in contrast, in- 
dependent of it as will become clear later. Insofar 
and with respect to the continuum limit, the vol- 
umes of the phase space elements are arbitrary 
and have no direct physical significance. Instead, 
physical significance is attributed to the mass 
rjjUJj (or the corresponding mass differences) of 
the phase space elements in //-space reflecting the 
properties of the basic counterparts LUj(r]jLUj) s in 
T-space. In this respect it is natural to incor- 
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porate the universal mass independence of the 
equations of motion of Violent Relaxation into 
the statistical mechanical picture by introduc- 
ing a constraint on the volumes Uj stating that 
the corresponding phase elements have constant 
mass, i.e. 

rjjujj = const. (6) 

We note that this procedure is in perfect agree- 
ment with the nature of the microscopic exclu- 
sion principle of Lynden-Bell (1967) which, for 
the case of only one phase space density, states 
that the mass of the phase elements is constant. 
The approach characterized by (||) may thus be 
considered to be just its generalization. How- 
ever, as a result and in contrast to Lynden-Bell 
(1967), where the phase space elements have dif- 
ferent mass but constant volume, the phase space 
elements considered here differ in volume while 
having constant mass. 

3. STATISTICAL MECHANICS OF VI- 
OLENT RELAXATION 

In this section we attempt to derive the coarse 
grained phase space distribution / of the final 
state of a system consisting of phase space ele- 
ments of J different phase space densities rjj sub- 
ject to Violent Relaxation. The densities rjj are 
obtained from the initial fine grained phase space 
density /. To determine the final state, we ap- 
ply the same maximum entropy principle as used 
by Lynden-Bell (1967). The state which max- 
imizes the entropy under the constraint of con- 
served total energy and mass is considered as the 
most probable final state attained by the system. 
As has been pointed out by several authors [e.g. 
Lynden-Bell (1967), Shu (1978), Madsen (1987), 
Shu (1987)], this final state is rather an ideal- 
ized final state since it implicitly relies on the 
assumption of complete mixing in phase space 
which eventually does not occur in real systems 
because of limited time scales. 

Let us start with the determination of the vol- 
ume tOj of a phase element of density m . Accord- 



ing to @, this volume is given by 

ujj = m/rij (7) 

where m is considered to be constant and is inter- 
preted as the individual mass of all phase space 
elements under consideration. In order to sim- 
plify the discussion but not introducing further 
restrictions, we now assume that the different 
phase space volumes Uj are all multiples of some 
smallest volume 

uj = mm{ujj} . (8) 

The volume Uj occupied by a phase space ele- 
ment of density rjj becomes 

Uj = gj u (9) 

where gj is the factor by which Uj is larger than 
the smallest volume uj, i.e. the microcell Uj is 
made up of gj microcells of elementary volume 
uj. 

The //-space is assumed to be divided into a 
large number of microcells of volume u each. 
At the macroscopic level, these microcells are 
grouped into macrocells containing a large num- 
ber v of microcells. The corresponding volume 
of the macrocell is vuj. Suppose that there is 
a microstate in which the ath macrocell contains 
n a j phase space elements of different densities rjj . 
Due to the collisionless nature of the interaction 
in the Violent Relaxation process described by 
the Vlasov equation (|T|) there is no cohabitation 
of microcells (Lynden-Bell 1967). 

We now calculate the number of ways of as- 
signing n a j microcells of volume uij and density 
rjj without cohabitation to the ath macrocell. 
Under the assumption that the phase space is 
resolved to the scale of the volume of the small- 
est microcell uj we find that the number of ways 
of assigning n a j microcells of volume uij without 
cohabitation to the ath macrocell is, according to 
(H), approximated by v\j{y — gjn a j)\. This num- 
ber includes cases where the microcells of volume 
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ujj assigned to the macrocell may not entirely be- 
long to it. However, these particular situations 
are rare since their relative occurrence is propor- 
tional to ~ l/sjv, where v 3> 1. In order to de- 
termine the total number of ways of assigning the 
J2j n aj phase space elements to the ath macrocell 
we assume the gj are ordered according to their 
values. Starting with the highest value assigned 
to gi or with the largest phase space element of 
volume ll>i, respectively, the total number of ways 
of assigning the J2j n aj phase space elements to 
the ath macrocell is 

,.i 

w(n aj ) 



{v - g\n a x)\ 

(u - gin a i)] 
({v - gin a x) - g2n a2 )l 

(" E; ■U J "„ J )\ 
iy - E/ 9j n < 



a j , 



Ei 9jn 



(10) 



0.3) 



By the ordering of the Oj and the exclusion 



principle expression (10) does not exclude every 
situation corresponding to partially overlapping 
phase space elements. Consider for instance a 
situation where a phase space element of volume 
oj j is assigned to a position between regions of 
the macrocell which have been already occupied. 
If the volume u>j is larger than the volume be- 
tween the occupied regions then a situation is re- 
covered in which phase space elements may par- 
tially overlap. However, as is the case for phase 
space elements overlapping the macrocell, situ- 
ations representing overlapping phase space ele- 
ments are also rare. This applies in particular to 
cases characterized by n a jUj <C ^,Vj. We will 
return to this question in a subsequent section. 
Here we assume that expression ( |i~0| ) yields a suf- 
ficiently good approximation to the total number 
of ways of assigning the J2j n aj phase space ele- 
ments to the ath macrocell. 

The total number of microstates W({n a j}) 
corresponding to a given set of occupation num- 



bers {n a j} is found by multiplying the expres- 
sions and taking into account the number of 
ways of splitting the total of N distinguishable 
elements into groups n a j. This yields 



JV! 



W{{n a3 }) = \{ , AJ 

j Ha n aj- a 



<->.) > 



n 



n 



.(11) 



Since the numbers w(n a j) are approximations, 
expression (|ll|) has also to be considered as an 
approximation . 

The macroscopic constraints the system is 
subject to are the J constraints related to con- 
servation of the total number Nj (or total mass 
Mj) of phase space elements of densities rjj 



E 'U-j""j = W N 3 = M 3 



The energy constraint reads 



(12) 



EEw«;(k + ^«) = £ (13) 



where the gravitational potential is defined by 

GrjjLOjTib 



<S>a = $(X Q 



E E 

j b=l,b^a 



(14) 



Note, that the volumes uj or, respectively, the as- 
sociated masses UjTjj enter the calculation through 
these macroscopic constraints. 

The most probable state is found by the stan- 
dard procedure of maximizing log W subject to 
the constraints of constant total energy and con- 
stant mass. Introducing the Lagrangian multi- 
pliers ctj and /3, respectively, related to the con- 
straints (|l^) and (|13|) the expression to be max- 
imized is 

ln[W({n aj })] -ay EE W n <y- 

3 a 

^EE%wk+k). (15) 
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The variation of dl3) leads to 



It is thus given as 



E I - ln [ n <y ] + - E 9jn a j]j Sn aj - 

n i E '/./-'/' ) ""./ ( 16 ) 

a 







and by taking into account the independence of 
the 5n aj 's, equation (|i~7|) splits into the J equa- 
tions 



In 



^ - Ej gjfi. 



a.l 



n 



"j 



a 



+ PVjUje a , (17) 



U = 1, • • • , J) 



where e a = ^v 2 + 3> a stands for the total energy 
per unit mass related to the ath microcell. The 
corresponding occupation numbers {n a j} are 

n a j = (y - 9jn a j) exp[-ajT]ju;j - Pr]ju;j€ a } . 

(18) 

By multiplying ( |l8| ) with gj , summing over j and 
subsequent algebraic manipulations we find 



(" - E 9j 



n 



aj , 



(19) 



J2j 9j expl-ajrijUj - j3r] j u: j e a \ + 1 



According to (p^), with fij = —ctj/f3 and with the 
constant ufflj = m from equation ([?]) the most 
probable occupation numbers {n a j} become 



71, 



a.) 



z/exp[-/3m(e a - fij)} 
Ej 9j exp[-/3m(e a - fij)] + 1 



(20) 



The coarse grained phase space distribution / is 
defined as the sum of the J phase space distribu- 
tions fj where 



/j(v,x) Fa fj(v a ,x a ) 



9j n ajVj 



(21) 



/(v,x) « /(v a ,x a ) 

= J2fj(Va,Xa) 



E 



9j n ajVj 



(22) 



Substituting for r\ a j from (|20j), the coarse grained 
phase space distribution / becomes finally 



/(v,x) 
E 



(23) 



9jVj exp[-/3m(e(v, x) - fij)} 
Ej ffj exp[-/3m(e(v,x) - //,)] + 1 



where e(v,x) = ^v 2 + ^(x) is the total energy 
per unit mass. According to the initial assump- 
tions, the state characterized by the occupation 
numbers ( p0[ ) or the coarse grained phase space 
distribution (|24|), respectively, is identified with 
the (idealized) final state attained by a system 
consisting of phase space elements of different 
densities rjj and which is subject to Violent Re- 
laxation. The Fermi-like functional form of the 
respective distributions ([2C|) and (24) originates 



from the exclusion principle which has been ap- 
plied in order to exclude situations corresponding 
to cohabitated microcells. 

One should keep in mind, that the distribu- 
tions ( |20"|) and (|4|), respectively, have been de- 
rived under the approximations made in evalu- 
ating the total number of microstates W({n a j}) 
corresponding to a given set of occupation num- 
bers {n a j}. However, for the two limiting cases 
(fj <C i]j , Vj) and (fj Fa r\j ,3j) the approxima- 
tion made in ( |l0|) is very good since for both 
situations the overlap of phase space elements of 
different volumes ojj does not occur. 

In finishing this section let us return to the 
question of the explicit determination of the mi- 
crocell and macrocell volumes coj and vu>, respec- 
tively. While the volumes Uj are absorbed in the 
common constant mass factor m = oojTij due to 
relation (0), equations ( 2C -p4]) illustrate by the 



canceling of v the arbitrariness of the volume vuj 
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of the macrocell. The only constraint thus re- 
maining is related to the applicability of statisti- 
cal methods, i.e. the constraint v ^> 1. 

4. DISCUSSION AND CONCLUSIONS 



In order to discuss the result ( p4j ) we compare 
the equilibrium distribution / with the coarse 
grained phase space distribution /lb which was 
originally derived by Lynden-Bell (1967) 



/lb(v,x) 



(24) 



E 



?? j exp[-/3 J (e(v,x) - ^-)] 



-f Ej exp[-Pj(e(y,x) - p,)] + 1 ' 

In the non-degenerate limit characterized by fj <C 
rjj , Vj, the coarse grained phase space distribu- 



tion (25) turns into a sum of Maxwellians of dif- 



ferent velocity dispersions <rj = /3- 

/lb(v,x) = ^2vj exp[-/3j(e(v,x) - /^)] (25) 

i 

where (3j cx In the non-degenerate limit, the 
final state of the process of Violent Relaxation 
is therefore identified with a thermalized state. 
This fact has been criticized by Madsen (1987) 
and Shu (1987). The velocity dispersions of the 
individual components are inversely proportional 
to their phase space densities and, as a conse- 
quence, the components have the same tempera- 
ture, a fact one would expect for a system sub- 
ject to a two-body relaxation process. This is not 
only counterintuitive, but it is also in disagree- 
ment with observations of astrophysical objects 
which are believed to have been subject to Vio- 
lent Relaxation. 

In contrast, the corresponding result derived 
form ( ]24| ) turns out to reflect the independence 
of the Violent Relaxation process on the phase 
space densities rjj of the different components 
involved. In the non-degenerate limit (fj <C 
r]j , Vj), the coarse grained phase space distribu- 
tion ( |24"|) does again become a sum of Maxwellians. 
But the Maxwellians are all characterized by 



one single velocity dispersion a 2 = ■ The 
coarse grained phase space distribution in the 
non-degenerate limit is 

/( v > x ) = E* ,? ?i ex P[~^ m ( e ( v ' x ) • ( 26 ) 

j 

We note, that according to Section ||] the approx- 
imation of W({n a j}) on which ( |26| ) is based is 



very good, because the condition n a juij <C u , Vj 
is equivalent to the condition fj <C r]j , Vj. In the 
non-degenerate limit, the final state of the Vio- 
lent Relaxation process defined by (|26| ) is the su- 
perposition of Maxwellians of a common velocity 
dispersion. This is in agreement with the com- 
mon velocity dispersion of the different compo- 
nents observed for instance in clusters of galaxies 
(e.g. Lubin & Bahcall 1993). The common ve- 
locity dispersion is equivalent to an equipartition 
of energy per mass. As a consequence, there is 
no mass-segregation present in systems subject 
to Violent Relaxation as has also been observed. 

Another limiting case for which the approx- 
imation made in evaluating the coarse grained 
phase space distribution (||) yields an exact re- 
sult is found for (fj rjj , 3j). This situation cor- 
responds to a dense coverage of the phase space 
by phase space elements with one specific density 
Vj- 

For the case of dense coverage of the phase 
space by phase space elements of different vol- 
umes u)j (i.e. fj~r]j,\/j) the situation becomes 
more complicated because the approximations 
which led to equation ([ll]) worsen. It is, how- 
ever, clear that the effective total number of mi- 
crostates corresponding to a given set of occupa- 
tion numbers {n a j} is less than W({n a j}) since, 
as discussed in Section S, the number W({n a j}) 
includes cases which correspond to partially over- 
lapping phase space elements. The question re- 
mains open of how the correction of this approx- 
imation would influence our result. An investi- 
gation of this point is left for subsequent work. 

In summary, we have performed a statistical 
mechanical investigation of Violent Relaxation 
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of phase space elements of different mass densi- 
ties. A well motivated constraint on the volume 
of the phase space elements has been introduced 
stating all phase space elements independent of 
their density to be of constant mass. As a conse- 
quence, the statistical mechanics accounts for the 
independence of the equations of motion in the 
process of Violent Relaxation on the mass density 
of the phase elements. As in the case of Lynden- 
Bell (1967), the final state reached by a system 
which is subject to Violent Relaxation is char- 
acterized by a Fermi-like coarse grained phase 
space distribution / reflecting the microscopic 
exclusion principle (Lynden-Bell 1967) that has 
been applied in order to account for the colli- 
sionless nature of Violent Relaxation. However, 
the functional form of this coarse grained phase 
space distribution is independent of the different 
phase space densities of the interacting compo- 
nents. In the non-degenerate limit, the coarse 
grained phase space distribution / becomes a 
superposition of Maxwellians, fi, all possessing 
the same velocity dispersion. This is in agree- 
ment with the observation of a common velocity 
dispersion of different components of astrophys- 
ical objects which are believed to have under- 
gone Violent Relaxation (e.g. the different mat- 
ter components in clusters of galaxies). Thus, the 
velocity dispersion problem which was present 
in the original treatment of Violent Relaxation 
(Lynden-Bell 1967) has been removed. 

As has been pointed out for a long time (e.g 
Lynden-Bell 1967, Shu 1978, Madsen 1987) and 
confirmed by N-body simulations (e.g Melott 
1983, McGlynn 1984, White 1996), Violent Re- 
laxation eventually fades before the final state 
(eq. |24}] ) is attained throughout a real astro- 
physical systems. However, expression ( |24| ) is 
likely to hold in the central regions where Violent 
Relaxation occurs most violently. In comparing 
the results obtained here with results from N- 
body simulations, one should probably take into 
account recently expressed scepticism concern- 
ing the reliability of N-body simulation. As has 



been shown (e.g. Kuhlman et al. 1996, Mad- 
sen 1996), there possibly exists a problem for 
high-resolution N-body simulations due to colli- 
sion and discreteness errors. When dealing with 
Violent Relaxation of nonbaryonic dark matter 
as, e.g., massive neutrinos (Kull et al. 1996), this 
problem becomes worse since massive neutrinos 
behave like an almost continuous medium. 

The assumptions and approximations on which 
the present investigation relies on are (1) that the 
process of Violent Relaxation can be described 
in the continuum limit, (2) that Violent Relax- 
ation of a multi-component system is essentially 
a collisionless process leading to a final equilib- 
rium state, (3) that the characteristics of the 
Vlasov equation in a statistical mechanical pic- 
ture of Violent Relaxation have not only to be 
represented by a microscopic exclusion principle 
but also by considering phase space elements of 
constant mass, and (4) that the approximations 
made in deriving the occupation numbers are 
valid at least in the non-degenerate case. As a 
consequence of condition (1), the present investi- 
gation in terms of phase space elements does not 
refer to its counterpart formulated in terms of 
particles (Shu 1978). The main difference com- 
pared to Lynden-Bell (1967) is the generalized 
constraint on the mass of the phase elements. 

The assumptions and approximations made 
here do not allow to consider the important case 
of incomplete Violent Relaxation. They also do 
not apply to the fact, that the final state reached 
by Violent Relaxation is, in a global picture, 
a transient state since two-body collisions will 
eventually start driving the system towards a 
thermalized state. In addition, the present inves- 
tigation does eventually not apply, as discussed, 
to all cases of dense coverages of the phase space 
by phase space elements. However, we believe 
that a more sophisticated treatment of the prob- 
lem which corrects for these approximations will 
not change our basic result. 
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